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, We consider the problem of the persistence of invariant curves for analytical fibered 

holomorphic transformations. We define a fibered rotation number associated to an invariant 
^ : curve. We show that an invariant curve with a prescribed fibered rotation number persists 

under small perturbations on the dynamics provided that the pair of rotation numbers verifies 
a Brjuno type arithmetical condition. 
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! 1 Introduction 

Skew-product transformations over irrational rotations have been vi^idely studied as a source of 
examples of interesting dynamics and for modehng fundamental phenomena in mathematical 
J> ' physics (c.f. quasi-periodic Schrodinger cocycle, see [1,14]). In the former direction, Furstenberg 

. [7] constructed the first example of a minimal non ergodic conservative diffeomorphism of the 

torus as a fibered circle diffeomorphism. Latter on, the basis for the study of these maps were 
CN ' stablished by M.Herman [8], who defined in particular the notion of the fibered rotation number 

(see also [11]). In the recent years, this theory was relaunched mainly by the works of G.Keller, 
\ J.Stark and T.Jager [9,10], who established (among other things) a Poincare-like classification 

relating the fibered rotation number to the existence of invariant graphs. 

In his doctoral thesis [21] (see also [22]), O.Sester studied hyperbolic fibered polynomi- 
als, successfully generalizing the classical notions of Julia set. Green's function, and the prin- 
^ \ cipal cardioid of the Mandelbrot set in the parameter space. Closely related works due to 

■ M. Jonsson [12, 13] are also important contributions to the subject. 

In this work we study fibered holomorphic dynamics. More precisely, given an irrational 
number a > 0, an open simply connected subset [/ of C, and a real number 5 > 0, we consider 
transformations of the form 

F : BsxU — > BsxC 

{e,z) ^ {e + a,fe{z)), 

where Bg denotes the strip |^ G C/Z | \Im{9)\ < 5}. We also assume that this transformation is 
a holomorphic function as a two variables complex function. In particular, the maps fe'.U — > C 
are holomorphic for every 9 £ Bg. We call F an analytic fibered holomorphic transformation, 
and we refer to it as a fht. 
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Note that a map F as above has neither fixed nor periodic point. The natural object that 
plays this role for the local dynamics is an invariant curve, that is, a holomorhic curve u : ^ U 
such that, for every 9 e Bg, 

F{9,u{e)) = {9 + a,u{9 + a)), (1) 

which is equivalent to fg(u{9)^ = u{9 + a). Indeed, the dynamics of F is organized around this 
"fibered fixed point", thus generalizing the role of a fixed point for the local dynamics of an 
holomorphic germ g : (C,0) (C,0) (see [16,18]). We define an appropriate notion of fibered 
rotation number, which appears to be quite useful for studing this dynamics. 

Fibered holomorphic dynamics may be thought of as a toy model for understanding the so 
called Melnikov's problem on the persistence of lower dimensional invariant tori for integrable 
Hamiltonian systems. Let us recall that, concerning this problem, celebrated results due to 
H.Eliasson [6] and J.Bourgain [3] establish the persistence of lower-dimensional invariant tori 
provided that the tangencial and normal frequencies verify a diophantine condition DC{t). 
However, it is conjectured that condition DC{t) is not optimal for this result. 

Inspired by this result and the corresponding conjecture, we consider the analogous problem 
for fhVs, namely the persistence of an invariant curve with prescribed fibered rotation number. 
The main theorem of this work states that invariant curves persist under small perturbations 
under a Brjuno type arithmetical condition Bi (which is more general than any condition of 
type DC{t)) between the rotation on the base a and the fibered rotation number (which in our 
setting corresponds to the tangencial and normal frequencies, respectively). 

Although //it's correspond to higher-dimensional systems, due to their skew-product struc- 
ture they have a one-dimensional flavor. In this direction, let us recall that the Brjuno condition 
is optimal for the linearization problem of both holomorphic germs and analytic circle diffeo- 
morphisms close to rotations (see [24,25]). In view of this, we conjecture that condition Bi is 
optimal for our problem. However, as in the classical case, the use of KAM techniques seems to 
be inapropriate to prove optimality. 

Let us finally point out that the problem of the persistence of invariant curves for C°° 
fibered holomorphic maps was settled by the author in [17], where it is proven that the optimal 
arithmetic condition in that setting is a DCi condition. 
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Doctoral, prepared at the Universite Paris-Sud XI. The author thanks Jean-Christophe Yoccoz 

for the suggestion of this subject, his advise and constant support to carry out this work. The 
author also thanks Jesus Zapata for a useful indication, as well as Jan Kiwi, Andres Navas 
and Godofredo lommi for several pertinent remarks and corrections. This work was partially 
supported by CONICYT-Chile and the French government. The article was written during a 
postdoctoral position of the author at Universidad Catolica de Chile and supported by this insti- 
tution and the PBCT-Conicyt via the Research Project ADI 17 on Low Dimensional Dynamics. 

2 Arithmetical conditions and linearized equation 

Trought this work, a > will be an irrational number (it will correspond to the rotation angle 
on the base for our fibered transformation). For a real number x denote by the distance to 
the nearest integer, that is, ||a:;|| = min^gz \x — p\- 
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2.1 Arithmetical conditions for a real number 

A complete treatement of what follows in this section may be found in [5,15,20]. For each iV G N 
we define the worst divisor of a up to the order N by 

Ta(N) = max ||na|| 

0<|n|<Ar" " 

Worst divisors may be completely described in terms of the denominators of the reduced of the 
continuous fraction of a. Indeed, the sequence go = ^,Qi,Q2, ... of these denominators has an 
important property, namely, for every n > 1, 

Qn+i =min{A: > qn | ||fcQ;|| < ||gna||}. 
Therefore, if A; G N is such that Qk < N < Qk+i, then Toc{N) = \\qk(^\\ ^- We define the set 

N>1 

If a belongs to B we say that a verifies the Brjuno (arithmetical) condition B (as defined by 
A. Brjuno, see [4]). According to H.Riissmann (see [19, §8]), this condition is equivalent to the 
convergence of the series 

^ log gn+l 
n>l ^" 

It is well known that B has full Lebesgue measure, but is a small set from the topological point of 
view (Baire's category). The following lemma contains a useful equivalence of Brjuno condition, 
which we will use through this work. 

Lemma 2.1 (Riissmann [19]) The number a belongs to B if and only if 

iogr«(2'^ 



2" 

n>0 



< oo. 



Proof . For every n > 2 one has 



^ i i + 1 ^ ^ i-\ V 

j=2" i=2" i=2" 



and therefore 



on + l _-) 
1 1 1 



E ,-2 



< > T7T < r-:;^n r< 



2"+i 2" 2"+i ^ i^ 2" - 1 2"+i - 1 2" ■ 

i=2" 



Since ra(-) is increasing, this implies that 



1 iogr„(2") ^ iogra(?;) ^ , r^(2"+^) 

i=2" 



which allows to conclude 



3 



2.2 Arithmetical conditions for a pair of real numbers 

Given a real number /3, we say that the pair (a, (3) is rational if there exists fc £ Z such that 
ka = (3 mod (1). Whe denote by Tj the set of non rational pairs. If (a, (5) belongs to Tj then, 
for a given G N, we define its worst divisor up to the order by 

Tcir{N)= max lino — ^. 

0<|n|<Ar" " 

This represents the quality of the approximations of j3 by multiples of a. Let Bi be the set 
defined by 

Bi = {(a,/3) G T| I a G B and ^"^^"/^""^ < oo}, 

n>l 

and let B" be the set of (3 such that (a, /?) lies in Bi. If a belongs to B then B" is a set of full 
Lebesgue measure which is meager (that is, equal to a countable union of empty interior closed 
sets). 

Lemma 2.2 Assuming that a lies in B, the pair (a, (3) belongs to Bi if and only if 

> < oo. 

n>0 

Proof . Analogous to the proof of Lemma 12.11 ■ 

2.3 Cohomological equation 

For more information on this classical subject see for instance [23]. Given an analytic function 
<j) : Bs ^ C, we look for an analytical solution -;/; : B5 ^ C to the cohomological equation 

i:{e + a)-^{d) = <i){6). (2) 

By integrating both sides of this equality, one readily checks that a necessary condition for the 
existence of a solution to this equation is that the mean value of (j) (with respect to the Lebesgue 
measure in C B^) is zero, that is, j^-^ (j){9)d6 = 0. A solution to the equation is not unique, 
but any two solutions differ by a constant. 

Proposition 2.3 If a belongs to B and the mean value of (p is zero, then there exists an analytic 
function ip: Bs ^ C which is a solution to ^ . Moreover, such a solution may be taken so that 

/^i i,{e)de = 

Although condition B is not optimal for this proposition, it is optimal for the problem of 
linearization of holomorphic germs in the neighborhood of indifferent irrational fixed points [24]. 



Remark 2.4 An easy homotopy type argument shows that the mean value of (j) on the circle T], 
is constant as a function of c £ (—6,6). If the value of this constant is non zero then we may 
envisage equation |^ replacing (p by (f) — Jji (j){6)d6 for some (any) c. 
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3 Definitions and normal forms 



In order to state properly the main theorem of this work we need to introduce some definitions. 

Let -F be a fhd and u an invariant curve for F. We consider the circle naturally embedded 
in Bs- In the remaining part of this work we will always suppose that f^i log \dzfe{u{9)^ \d9 = 0. 
(Notice that, since F is injective, the differential dzfe is always non zero.) This condition says 
that the curve is neither attracting nor repulsive at the infinitesimal level. We say that the curve 
is indifferent. We also suppose that the application 9 i-^ dzfe{u{9)) is homotopic in C \ {0} to 
a constant. 

Definition 3.1 We define the fibered rotation number by 

Qtriu) = j^\ogdzfe{u{e))d9. 

This number represents the average rotation speed of the dynamics around the invariant curve. 
Notice that the log above is well defined mod 2iTi, and the number Qtr{u) is well defined mod 1. 

Remark 3.2 Notice that, since all of the circles = |0 G i?^ | Im{9) = c} are homotopic 
between them and F is holomorphic, one has 

Qtrin) = j^\ogdzf{9M9))d9 

for every circle T].. Thus, the indifferent nature of the curve (resp. the fibered rotation number) 
may be detected (resp. computed) on any circle TJ, for \c\ < 5. 

Along this work we will often deal with holomorphic changes of coordinates H defined from 
a tubular neighborhood of the zero section {z = 0} to a tubular neighborhood of the invariant 
curve, say 

{9,z)^{9M{z)). 

The functions kg will be biholomorphic transformations between two topological discs sending 
the origin to the curve, that is, /i6)(0) = u{9). Moreover, the application 9 ^ /i6»(0) will be 
homotopic to a constant in C \ {0} (in other words, it will have zero topological degree). There- 
fore, when conjugating our transformation F by such a H we will get a new transformation 
F = H^^ o F o H having the zero section u = {z = 0}bs as an indifferent invariant curve 
and having the same fibered rotation number, that is, Qtr{u) = Qtr{u). In particular, the indif- 
ferent nature of the curve and the fibered rotation number are invariant by this type of conjugacy. 

Assume F : Bs x U Bs x C is a fhd, a is a Brjuno number, and u is an indifferent invariant 
curve for F. For /? = Qtr{u) we will solve the (cohomological) equation 

v{9) - v{9 + a) = 2mP - log dzfe{u{9)) 

on the strip Bs (see Proposition 12. Then letting ui = exp{v) and H{9, z) = (9, u{9) + ui{9)z^ , 
and performing the corresponding change of coordinates, one gets the following normal form for 
F: 

Nf{z, 9) = H-^oFo H{9, z) = {9 + a, e^^'l^z + p{9, z)) . 
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Here p is an analytical function vanishing up to the second order at z = 0. This function p, 
and so the flid Np, is defined on the strip Bg for the 6 variable, and p{6, •) is defined on a 
non-constant open set Uq = H''^{{9} x U). 

Assume now that F is a fhd that may be putted in a normal form 
H-^oFo H-\e, z) = {e + a, e^""'^ z + p{e, z)) 
by a change of coordinates of the form 

H{e,z) = {e,uo{e) + ui{e)z), 

where uq ■ Bs' ^ U, ui : By C are analytical functions for some 5' G (0,(5], and ui has zero 
topological degree. In this case one can readily check that uq is an indifferent invariant curve, 
and its fibered rotation number equals q{uq) = /3. 

Therefore, the existence of an indifferent invariant curve with transversal rotation number (3 
is equivalent, under the Brjuno condition over a, to the existence of an analytical fibered affine 
change of coordinates that puts F in an apropriate normal form. 



4 Statement of the problem and result 

Let F be a fhd with an indifferent invariant curve having a given fibered rotation number, say 
/? G M. As we claimed in the previous Section, if a G IB then F can be written in a neighborhood 
of the curve in a normal form {9 + a, e^'^'^z + p{6,z)), where p is an analytic function. The 
function p{6, •) vanishes up to the order 2 at 2; = 0, and is convergent for \z\ < r for some positive 
radius r. A small perturbation of such a transformation is defined as being a fhd of the form 

{6, z)^{e + a, poie) + (e^-^^ + + z)) , 

where po,Pi are small analytic functions, where p is an analytic function whose em size is 
comparable to that of p, and where p{6, •) vanishes up to the order 2 at z = and is convergent 
for l^l < r. 

Notation for sizes. Every function that will appear in this work is a bounded holomorphic 
function with respect to its variables (defined in some complex open set). Let f -.U dC^ ^ C 
be a bounded holomorphic function from the open set U C C*^, k G {1,2,3}, to C. We define 
the norm ||/||;7 of / as 

||/||t/ = sup|/(z)|. 

In general, U can be thought as a product B^ x x , where D^, _D* are discs in the complex 
plane. We then may define the oscillation of / with respect to the z variable as 

^^^if) B, DZ D* = ™P Ifi^^^i^t) - f{e,Z2,t)\. 

' {e,t)£BsxD* 

Note that if z = belongs to and / vanishes at this point, then its oscillation is a bound for 
its norm, that is, 

\\f\\Bs,Dz,Dt < OSc{f)j^^ j^^ j^^ < 2\\f\\B,^DZ,D^- 

We will currently omit the corresponding domains of definition when this does not lead to 
confussion. 
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4.1 1— parameter families 



In general KAM results, a perturbation over an elliptical dynamics (rotations of the circle, com- 
pletely integrable Hamiltonians, etc.) leads to a perturbation of the frequencies which control 
the dynamics. Thus, in the perturbed situation, we can not expect to retrieve the same dy- 
namical properties than in the unperturbed elliptic situation. In order to retrieve these original 
properties, one introduces a 1— parameter correction, where the parameter is a vector value of 
the same nature (and dimension) as the related frequencies. In this way, one usually shows that 
some dynamical properties (linearization, existence of invariant curves with a given rotation 
number, etc.) are persistent in codimension 1 (see for instance [2]). In our work, we will only 
perturb the holomorphic part of the transformations. Therefore, our frequency corresponds to 
the fibered rotation number, that is, a complex number. 

Let F be a fhd and u an indifferent invariant curve, with Qtr{u) = /? G M. Let S be an open 
set in C. A transversal small perturbation of F is a 1-complex parameter family {Ft}tgs of //id's 
(i.e., a complex curve in the space of //id's) verifying the following properties: 

- every element Ft is a small perturbation of F, 

- the fibered rotation number (even if an invariant curve does not exist) changes along the 
complex curve {-Ffjtgs (see Theorem 14. II for more explanations). 

We say that the curve u is persistent if for any transversal small perturbation {Ft}t(^Y, 
there exists a parameter € S such that Ff has an invariant indifferent curve u* satisfying 
Qtr{u*) = Roughly speaking, our result says that, except for a complex correction, any small 
perturbation of F has an indifferent invariant curve having /? as its fibered rotation number. 
However, an arithmetical condition over the pair (a, (5) is required. 

We say that a 1-parameter family Ft of //id's is analytic if the application 



is holomorphic. We consider the usual development for the elements of this family, namely. 



To fix ideas, we will suppose that pt{0, •) is convergent in the unit disc B. With this notation 
we can finally state our main result. 

Theorem 4.1 For every pair {a, (3) verifying the Brjuno condition Mi, and for every positive 
constants L > 1,M,T and 5, there exists a real number e*{L,M,T,5,{a, (3)) > and a positive 
constant Kji{L,M,T,5,{a, f3)) such that the following holds: if e belongs to (0,e*] and for an 
analytic family {i*t}tgs there exists a disc D{tQ, Kre) C S such that 



T.X Bsx D{0,r) 
{t,9,z) 



C 



Ft{e,z) 



Ft{e,z) = [e + a,po,t(.o) + {e^^"' + pi{e))z + pt{e,z)). 



(3) 




(4) 



and such that, for every t in D{tQ, Kpte) 
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then there exists a parameter t in D{tQ, K^e) such that Fj has an indifferent invariant curve u 
which is analytic on the strip Bs. Moreover, the fibered rotation number of this curve equals (3. 

2 

Furthermore, the size of u goes to when e goes to 0. 

Note that the transversahty condition is given by (jlj). In [17] we show an analogous result in the 
case of C°° transformations. Indeed, we show that a diophantine condition over the pair (a, (3) is 
sufficient in order to get persistence of the invariant curve. We also show that this arithmetical 
condition is optimal. For a non-diophantine pair (a, /?) we present a smooth fhd having a non- 
persistent indifferent invariant curve with (3 as fibered rotation number. We conjecture that in 
the analytic setting, the Brjuno Bi arithmetical condition is optimal in this sense. 

Plan of the proof. The proof of Theorem 14.11 will be obtained as an application of the 
successive conjugacy method (Newton's algorithm). We will find an analytic fibered affine 
coordinates change that will put one of the //id's Ft in the family in the adequate normal form. 
As usual in this type of proofs, we must keep the control on several quntities at the same time 
and perform many different operations. Thus, the whole proof is divided in various parts that 
we describe now. In Section we make a review of some simple but useful results in complex 
analysis that will be used in the demonstration of Theorem 14.11 Some conventions about the 
names of the constants arising in the work are done in that Section. In Section [H] we show two 
lemmas which represents the core of the proof. These lemmas define the convergence speeds of 
the iterative process. Moreover, they give the size of the width strip loses that we must introduce 
at each stage. These lemmas establish the relations between the arithmetical properties of the 
rotations numbers and the characteristics of the iterative process to be applied, and therefore, 
between arithmetics and dynamics. In Section [7] we introduce the actual iterative process. In 
17.11 we define a constant which determines the values given by the lemmas in Section [6l In 
17.21 we make a detailed description of the algorithm. We introduce and describe the different 
operations used at each stage. We also list the order in which these operations will be applied. 
In l7.3l we make some preliminary estimates and we assume two additional hypothesis in order to 
reach the adequate situation which enables us to start with the iterative process. In l7.4|7.5t 17.61 
we perform in a detailed way all the operations and estimates of a single stage of the algorithm. 
We show indeed that these operations allow us to get the good estimates in order to iterate the 
process and thus, to approximate the normal form better and better. This is the more delicate 
part of the proof, since we must control many quantities that arise as the result of the different 
operations. In 17.71 we show that our algorithm converges, and thus, it gives the desired normal 
form for one of the elements Ft in the family. In l7.8l we exhibit a modification of the first stage of 
the algorithm, in order to obtain more precise results when e is very small. In [8] we remove the 
additional hypothesis introduced in Section [7.31 We show that we can get the desired situation 
as a consequence of the original hypothesis of Theorem 14.11 by means of a previous preparative 
procedure. In [U] we prove a parametrized version of the main theorem. We expect that this 
version of the theorem should be useful in many other applications. 

5 Some results on analysis 
5.1 Constants 

In this work, several positive constants with different origin will appear. Some of them have a 
fixed numerical value. These constants are called universal and are denoted by Cj, where the 
index represents the order of apparition of the constant in the work. If we denote a constant 



8 



by Cj then this means that this constant is universal. There exists constants which make part 
of the statement of the theorem, these constants are L,M,T,6, In a third group, there 

are constants whose values depend on the values of the constants L,M,T,S,{a,/3), and these 
arc denoted by Kj. The index represents the order of apparition or the feature of the constant. 
The special notation H is reserved to a particular constant belonging to this last group. At some 
extent, this constant determines completely the proof. 



5.2 Estimations for truncate analytical functions 

Let / be an analytical function defined on the strip Bg, and let TV be a natural number. We 
define respectively the truncation up to the order N and the rest up to the order AT of / by 

f\N{e) = Yl 

\n\<N 
|n|>Ar 

where /(n) is the ra— th coefficient of the Fourier series of /. These two functions are analytical 
in the strip Bg. If S' e (0, S) and S — S' is small then there exists a constant Ci such that 

|/(n)| < ||/||B,e-2-l"l^ (5) 

\\f\w\K < CiWfWs, . (6) 

The last inequality says that a small loss in the width of the definition strip gives an estimate 
on the norm of the rest of / (Cauchy estimate). We recall here the following simple fact. 

AfRrmation 5.1 /// is an analytical function whose differential dgf is hounded in Bg then 

f{9)d9 + \\def\\B,. (7) 



5.3 Estimates for the solutions of some truncate cohomological equations 

We consider the following cohomological equations on gi and g2 

9i{0) - 9i{e + a) = pi{e) (8) 
e2-/352W-52(e + «) = P2(.0), (9) 

where a, (3 are real numbers, pi,P2 are trigonometric polynomials of orders smaller than N + 1, 
and jjipi{9)d6 = 0. Using the Fourier coefficients method we obtain the expressions for the 
solutions (as explained in the Introduction, for the first equation we just retain the solution with 
zero mean): 



0<\n\<N 



g27rma 



92 



(0) = E 



\n\<N 



P2(n)e2'^^"^ 

Q2iTif3 g— 27rma ' 
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Let /i,/2 be analytic functions defined on Bs, with Jji fi = 0. For j... we put pj = fj\^, for 
j e {1,2}. Denoting (pi{N) = TaiN),if2iN) = r„^^(iV) (see Section [2]), we get a constant C2 
sucli that 

||<7,I|b, < C2iV||/,||B,<^,(iV) 

for j € {1, 2}. Note that the factor C2N is far from being optimal, but it will be enough for our 
purposes. 



6 Two technical lemmas 

In this Section we show two technical lemmas which will provide the definition of two sequences 
of positive numbers {ln},{'Wn} converging to in a very particular way. These sequences will 
be used as a measure of the rate of convergence of our procedure. These two lemmas will 
also provide two summable sequences of positive numbers {dj^}, {d^}. These sequences will 
measure loses in the width of the analyticity strip in order to get some Cauchy estimates. This 
Section is the most important part from the point of view of the relations between dynamics 
and arithmetic. Let H > 2 be a real parameter and {^n}n>o, {-Bn}n>o be the non-decreasing 
sequences defined by An = Ta,/3{2^) and i?„ = Ta{2^). Under the Brjuno Bi condition these 
sequences verify 

log An ^ ^ log Bn ^ 

> < CXD, > < 00. 



Lemma 6.1 For 

the following holds: 

1. lim„^oo^n = 0. 

2. j:^m^Bnln < I- 



e 



-n2 



16i<AnBn2'^ 



(10) 



3. ^^2 Anln ^ \Q • 

4- The relation 



_ Mne 



(11) 



"n 



defines a sequence {d^}n>o of real positive numbers such that 

^4 <Ki(l+logH) (12) 

n 

for a constant Ki depending only on {A„},{i?„}. 
Proof. The proof of the statements 1., 2., 3. are straightforward. For 4- we put 

an = -log(H2'^/„) 



log An + log Bn + log 16 + n 



2 
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and we define a sequence {dl^}n>o by 
We obtain that 

27r2"d^ = <T„+i -a„ + log2H. 

As an is increasing, the numbers are positive. Since the series is convergent, the series 

^ dn is also convergent, and its value is bounded by Ki(l + log^^), where Ki only depends on 

{An},{Bn}. 

We construct now a sequence {a;^}„>o in such a way that the perturbed numbers = d^+x^ 
verify point 4.. For this, we define as being the solution of the equation 



d^ + 



1. (13) 



It is easy to see that (fT3l) has a unique solution x^. Moreover, we have x„ < s^, where s„ is the 
solution to 

This equality easily implies that the series ^ s„ converges, thus concluding the proof ■ 
Lemma 6.2 For 

^-=^^' 

where the sequence {ln}n>o is defined as in Lemraa \m\ the following holds: 
1. Hu;„(2"yl„)2 < ^. 

3. The relation 

— = <¥. (i^> 

defines a sequence {d^}n>o of positive numbers such that 

<i^2(l + logH), (16) 

n 

/or a constant K2 depending only on {An},{Bn}. 

Proof . The proof of the statements 1., 2. are straightforward. Let us show 3.. As in the proof 
of Lemma Em we define a sequence {d^}n>o by the relation 

that is, 

4 = ^{2^ + 1 + log(^^+ii?n+i) - \og{AlBn) + log32H}. 

Therefore, we have a positive and summable sequence. As in Lemma 16. H we can find a real 
sequence {x^}„>o such that the perturbed numbers d^ = d^ + verify 3. ■ 
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7 An iterative process 

7.1 The choice of K 

The sequences defined in the previous Section depend on the value of the constant K. In order 
to use these sequences in our process we need to choose the value being larger than many 
quantities which depend only on the values of constants of type K (cf. Section [5. This choice 
can be made a priori since the type K constants depend only on the values of the constants 
L,M,T,S, (a, /?) that appear in the statement of Theorem 14. li However, along the process we 
will write explicitly the bounds assumed for b^. Fixing the value of we obtain two positive 
sequences {d^}n>o, {d^}n>o- We will use these sequences as loses in the width of the analyticity 
strip. These sequences verify 

^max((i°,4) < (Ki + K2)(l + log^^) < oo- 

n>0 

We choose a natural number n* such that 

Yl ™ax(d°,4) < ^, 

n>n* 

where 6 is the initial width of the strip for the 9 variable. The proof of Theorem 14.11 will be 
made using an algorithm divided in stages. The number n* represents the starting stage for 
our algorithm. At each stage the functions will be analytic on a strip (for the 9 variable) whose 
width equals 

<5"* = - (17) 
4 

^n+i ^ 5"-max(d°,di). (18) 
The choice of n* implies that for every n > n* one has J" > |. 

7.2 Description of the algorithm 

The proof of Theorem 14.11 will be obtained by showing that there exists an analytic conjugacy 
of the form 

h{9,z) = {9,uo{9)+ui{9)z) 

that puts one of the Ft into the adequate normal form, with ui having zero degree. (As noted 
in Sectional this implies the desired result.) After conjugacy one obtains 

h-^oFoh = (^9 + a,^^j^{po,t{9) + Xuo{9)-uo{9 + a)} 
+ T^^{pi,mMO) + Pt(.9,uo{9))} 

+ z^^{X + p^,t{9)+dM9,nom} 

+ H^ika) {Pti(^^ MO) + uii9)z) - pt{9, u,{9)) 

- zui{9)d,pt{9,UQ{9))}). 

We will perform a process in order to eliminate the terms /Oq, Pi (this implies the desired normal 
form). This process is divided in many stages. The goal of each stage is to decrease the sizes of 
Pq and pi, in such a way that at the limit these sizes are 0. However, the rate of convergence is 
very slow and deeply depends on the arithmetical conditions for the pair (a, /9). 
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At each stage we will apply four types of operations. The first two ones are KAM type (or 
Newton algorithm type) operations. That is, we solve the related equations but in a simplified 
way (we only keep terms which turns these equations into cohomological equations). These 
operations provide the functions uq,ui as solutions to cohomological equations, which give the 
corresponding coordinate change related to the stage. Obviously, this coordinate change does 
not conjugate our transformation into the desired normal form, since we only solve simplified 
equations. We must deal with some rests, whose sizes are controlled via the remaining two 
operations. The third operation consists on introducing a loss in the width of the analyticity 
strip for the 9 variable. The fourth operation consists on restricting the parameter space. 

We describe now the first two KAM type operations. These operations consist on solving a 
truncate cohomological equation (cf. Section [53]) : if we are in the stage n we solve the equations 
with truncation up to the order 2". 

7.2.1 The four operations 

Case ui = I (solving the uq equation). After conjugacy by a map h for which ui = 1 one 
obtains 

h-^oFoh = (^e + a,po,m + Xuo{e)-uo{9 + a) + pi,t{e)uo{e) + pt{e,uo{e)) 

+ z{X + pi,t{0) + d,pt{9,uo{9))} 
+ {pt{9, uo{9) +z)- pt{9, uo{9)) - zd,pt{9, uo{9))} 

Choosing uq so that 



Po,m + Xno{9) - uo{9 + a) = 0, (19) 

2" 



the new dynamics is given by 



2" 



PoT = + Pt{-,uo) + po,t 

PiT = Pi,t + d.Pt{;uo) 
pT'^i-^z) = pt{-,uo + z) - pt{-,uo) - zdzPt{-,uo). 



A crucial estimate will be 



Inoll <C22'^r„,^(2")||/5o,t||. (20) 



Having in mind the sizes of po.t and pi^t, using this inequality we will see that the size of 
puew _ decreases ("quadratic" effect). The sizes of the remaining involved functions 

remain almost unchanged. 

Case Uq = (solving the ui equation). After conjugacy by a map h for which uq = one 
obtains 



h-'oFoh = [e + a,£^ + z^{X + p,AO)} 

Pt(e,uiie)z) \ 

ui{e+a) )■ 

It would be natural to choose ui so that 

^i^^{A + pi,(m = A (21) 
ui[9 + a) 



13 



for an adequate constant A (given by the zero mean condition in the cohomological equation). 
However, due to technical reasons, we will only deal with the corresponding 2n truncated coho- 
molical equation The resulting dynamics is then given by 

PoA-) 
ui{- + a) 

A - A + Rest{0,t) 
pt{-,ui{-)z) 
ui{- + a) 

where the expression of Rest{9, t) will be made precise in the next Sections. For ui we will get 
an estimate of the form 

\\ui-l\\<i<2^T^{2^)\\depi,t\\. (22) 

The bound for the sizes of pQf" and pf^^ remains essentially unchanged. The effect of solving 
an equation as the truncated version of ()2ip may be interpreted as an attempt for turning pi 
independent of 6. Indeed, this operation allows us to control the size of dgpi^f However, to 
accomplish this task we must control the size of Rest{9, z). 

Loss in the width of the strip. The third operation consists on applying Cauchy estimates 
on the sizes of the functions. For this, we take the sup on a smaller strip. We call this operation 
a loss on the analyticity strip. By this way, we control the resulting rest from the two precedent 
operations. At the stage n we will lose max((i^, d^). The new width of the strip will be as in 
p8|) . This operation provides a control on the size of Pq^" and dgp^^^. 

Reduction of the parameter space. Due to Claim [SHI in order to control the size of p^^^^ 
we need to control the size of dep^f" and the complex number Jji p^^^ {9)d6 . The latter complex 
number may be changed with the parameter t due to the transversality condition on the family 
{Ft}te'E- Thus, this operation consists on localizing a disc in the parameter space where this 
complex number is small enough. We pick a disc centered at a simple zero of Jrj-i p^^^ {0)d9 , and 
with a very small radius. 

7.2.2 Description of a stage 

This Section only presents an informal description, and the complete details will be postponed 
until Section EH At the beginning of the n— th stage we have the functions /9g^,p"j,p", and 
the estimates 

I|P0,J < (23) 
IIPmII < i^3/n (24) 

for a constant ii'3. The sequences Z„, Wn are those defined in Section [6l We also have a parameter 
space disc D{tn,Pn)- The goal at the end of the stage is to achieve the corresponding necessary 
estimates to begin with the n + 1-th stage. Each stage splits into 4 parts: 

First part. In this part we apply many times the "solving the uq equation" operation. 
Hence, this part is composed by many steps. We introduce a supplementary index indicating 
the corresponding step. We will obtain the functions p^'l'^^ , p^'l'^^ , p^'^~^^ as the result of the 
step i. At each step the sizes of the functions pi^tiPt remain essentially unchanged. Function 



jnew 

Po,t 



„new 
Pl,t 



new 
Pt 
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Pq'I is the sum of a quadratic term whose size decrease in a half at each step, plus the series of 
rests having orders greater than 2" (Fourier series orders). The first part ends when the size of 
the quadratic term is smaller than '^"+^ . 

Second part. This part consists in applying one time the "solving the ui equation" operation. 
This enables to relate the control on the size of d^pi^t to a control of a rest having order greater 
than 2". 

Third part. This part consists in introducing a loss in the width of the strip. The first two 
parts provide some rests arising essentially from truncate functions whose sizes are smaller than 
Wn,ln, respectively. Inequality the relations (dH), and a loss on the strip width of an 
amount of max((i^, d^), provide the estimates UPo^^H < Wn+i, WdePi^^W < In+i- 

Fourth part. This part consists in applying the operation which reduces the parameter 
space. We get the bound In+i for the complex number Jji Pi~^^- Hence, we obtain all the 
necessary estimates to start with the n + 1-th stage. 

Therefore, the proof of Theorem l4. II depends on the possibility of showing that the hypotheses 
on the family {Ft}t(^s allows, on the one hand, to obtain good estimates to start the first stage 
(stage n*), and on the other hand, to complete our plan for any stage. We will see in Section [7.71 
that the properties ensured by the Lemmas 16.11 and 16.21 for the sequences ln-,Wn allow to show, 
for only one parameter t* in S, the convergence of our process. In other words, the successive 
composition of coordinates changes provided by each stage converges to an analytical coordinates 
change which puts the fibered holomorphic transformation Ff in an adequate normal form. 

7.3 Initializing the algorithm 

At the beginning of our process we will deal with functions defined on the following domains: 
the unitary disc B = Z)(0, 1) C C for the z variable; the strip of width 5 for the 6 variable; 
and the disc I?(0, K'e) C C for the parameter t, where the constant K' will be explicit very soon 
(see equation (I29p ). Before starting the iterative process we need to estimate the derivatives 
with respect to 6 of some functions. Allowing a loss of 5/4 in the width of the strip, classical 
Cauchy estimates give 

24M 

Wded.PtW < ^=--N, 
d 

where the last equality provides the definition of the constant N. We may already exhibit the 
minimal value for the starting sizes, namely, 

e-^. (25) 

These sizes only depend on 6 and n* , which in their turn only depend on L, M, T, 5, {a, (3). We 
will assume in this Section that 

WpoW < Wn*. (26) 
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Notice that this hypothesis is stronger than those required in the statement of Theorem 14.11 
However, we will see that we can reduce the general case to this one by a preliminary process 
whose description we delay up to Section 18.11 We consider the functions 



n*,0 



PO,t 



n*,0 



Pl,t 



n*,0 
Pt 



Pt 



(27) 



defined on the following domain: the strip of width 5" = 3(^/4 for 9, and the disc D{0,R^ 
of radius i?" '"^ = 1 for z. The parameter space is the disc D{tn* ,Pn*), centered at t„* = 0. We 
will assume that tn* is a simple zero of J^i pl^ . (This hypothesis will be removed at Section 
The radius of the parameter space disc is 



Pn* 



100/„* 



2400e 



(28) 



Hence we pick 



K' 



2400 
~5~' 



(29) 



Since H is a large number, we have that p„* < 1. Up so far, we have the following estimates to 
start with the stage n* : 



115.^/11 


< 




(30) 


\\po7\\ 


< 


Wn* 


(31) 


\\dzPt II 


< 


M 


(32) 


osc{d0dzPt '°) 


< 


2N 


(33) 


osc{dtd^p'} '°) 


< 


2T 


(34) 




< 


L-i 
1000' 


(35) 



Notice that inequality (j35|) and the existence of the complex number Aq are also additional 
hypothesis which we will assume for a while. (These hypothesis will be removed in Section [HI) 
Finally, we assume that L > \Aq\ > L^^. 



Remark 7.1 (On the starting stage and size of e) We will see in Section 7.8 that a subtle 
modification in the first stage yields the last claim of Theorem \4-ll which concerns the sizes of 
u and £. 



7.4 Realization of the iterative scheme 

To control the many functions arising during the realization of the iterative scheme, we introduce 
some real sequences. (The justification of some claims on these sequences made below will be 
deferred until they become transparent from the point of view of the iterative process.) At the 
begining of the stage n we have the family {Ft} written in the form 

F:^\e, z) = {e + a, plf{e) + zpif{e) + \z + p^'^e, z)) , 

where p"'" vanishes up to the order 2 at z = 0. The definition domain for ^p"'" is the disc 
D(0,i?"''^) for the z variable, where the sequence {-R"'*}n>n*,jgNu{oo} verifies 
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and the strip B^n of width (5" defined by (llSp for the 9 variable. The parameter space is the 
disc D{tn,Pn) which is centered at a simple zero tn of Jji Pif{9)d9, and whose radius equals 

Pn = WOln- (37) 

In this Section we will deal only with the estimates which do not concern the parameter t, and 
we defer to Section [7.51 the estimates involving this parameter. At the beginning of the stage we 
have the estimates 

(38) 

(39) 
(40) 
(41) 

where the real sequences {Mn}n>n* and {-/Vn,i}n>„*,igNu{oo} are bounded from above by some 
constants Km and Kjy respectively. 

Using Claim [STTl inequality (j4ip . and an estimate on the norm || J^i dtPif — Ao|| (see in- 
equality (pll) in Section 17. 5p , we obtain 



WPofW 


< 


Wn 


\\dbt''\\ 


< 


Mn 


osc{dedzpl' ) 


< 


Nnfl 


\\depl^!\\ 


< 


1 



IIPiJII < K^ln (42) 

for some constant K^, > 0. 

7.4.1 First part: the uq equation 

Suppose that at the first part of the stage n we have the functions Pq'I, p^'l, and p"'*, where the 
last one is defined on the disc D{0,R"''^) for the z variable. We will describe the i-th step. In 
this step we search for a coordinate change of the form h\^^{9,z) = {9,Uq'1{9) + z) by solving 
the truncate equation 

^^0,1 (^) - ""o,* (6' + a) = -Poll 1 2" • (43) 
By the results of Section 15.31 the solution Uq'I satisfies 

H:i\\ < C22"r„,^(2«)||;.^;;|| (44) 

Wdeu^o'l < C24"r„,^(2")||p^;;||. (45) 
Notice that the resulting new po is 

n,i+l n,j| . n,i n,i , nA / n,i\ / An\ 

Po,t = Po,t\w + Pi,t'^o,t + Pt {-^UQi), (46) 
which due to the definition of PQ'l\2r^ satisfies 

We also get a new pi, 

nA+l n.i , nA/ nA\ /ac>\ 

Pi,t = Pi,t +dzPt {-jUq ), (48) 

as well as a new p, 

nA+l/ \ riA/ nA , \ nA / nA\ n nA/ nA\ t ac\\ 

Pt {-^z) = p^' {■,Uq' + z) - pt' {■,Uq') - zd^pt' {■,Uq'). (49) 
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The function _r'^((9, •) is defined on Z)(0, r"''+\6i)), with 

r^''+\e)=R^'^-Ki{e)\. (50) 



The sequence {i?"'*} verifies 



< r'^'^+^e) (51) 



for every 6 in Bs^- From now on the domain for z is D(0, R^'^^^). We recall that equation (j43p 
is solved without lossing width on the analyticity strip. In fact, we will introduce a loss of width 
of size (Iq only at the moment when we will need to estimate the rests j l^^r- This will be made 
only once per stage, and at the third part of the each stage. 

Estimates at the first part of the stage. We define 

Vnfi = Po,t (52) 

Vn,i+1 = Pl,tUo,t+ Pt i^^UQ't). (53) 

Since (|17|) yields 

n,i+l\ t I /rA\ 

P0,t \2" = Vn,i+l\2" , (54) 

each time that we solve equation (I43p we actually obtain a better bound on r/^ This allows 
to improve dH]) and ([I^ to 

Wu^W < C22"r„,^(2")||<J| (55) 
Wdou^oiW < C24"r„,^(2")||<,||. (56) 
We will show by induction on the index i > of the corresponding step that, for some constant 

(57) 
(58) 

(59) 

(60) 

(61) 

During the proof we will also provide an inductive definition for the sequences {i?"''*}j>o and 
{]^n,i}i>o- Notice that inequalities ([5^ . (fSHj) and ([SU]) give 

||^".;|| < ^MlL (62) 
II o,tii - ^2"2* ^ ^ 

Wdeu^ltW < (63) 



II o2 n,iii 

\\dzPt II 


< 


Mn 




< 


Nn,i 


II ^n,i II 


< 


Wn 
2» 


\\9ePi;l\\ 


< 


2/ 


\\Pi';\\ 


< 


Kiln 



The case i = for the induction is furnished by the hypothesis (|38l) . . . . . (I42|l . We now assume 
that the desired relations hold for every j G {0, . . . , i} and that we have at our disposal the 
constants R^'^ and Nn^j for every j G {0, . . . , i}. We let 

^ i?".'-C22"r„,^(2")t^„(i) (64) 
^n,oo ^ i?"'0 _ 2C22"r„,^(2")w;„. (65) 



18 



If is larger than C2 then inequality (I5ip follows from ()50p and (I62p . Inequality (I6ip follows 
from (08]), (021), dSZD and dMD, since 



l/'u II < 



j=0 



< K3ln + MnY,\K,i\\ 

j=0 

< Kiln. 

Notice that the constant Ki^ only depends on C2,i^3 and Km- From relation (09]) we obtain 

RV'^+^ii = Rvr(-,^+<'*)ii<M., 

which implies (|57p . Differentiating (09]) we get 



eOzPt = ogd^pt ' {■,z + Uq'^) + d^pt ' {■,z + u^'^JdeUf^ 



n,i\ . o2 / 



sdzPt «o t ) - ^zPt (•> «o < )^e^o 



By estimating the corresponding terms we can bound the oscillation 



osc{dedzp'l''^') < iV„,. + 2i^Mi^24"r„,^(2-)|f. 

The above inequality allows us to give the definition of the sequence {Nn,i\i>Q'- letting A^^^'^o 
2N and assuming that N^^i is already defined, we put 



We also define Nn,oa as the supremum of the Nn/s, that is, 

Nn,oo = Nn,0 + 4KMC24"r,,^(2")t«„. 

It is not hard to check that 

osc{dedzP^''^^) < NnA+i- 
Relation (08]) allows us to (compute and) estimate dep^'^^^ by 



(66) 



(67) 



dn.i+l I 



dep'lf + Y.MdzPt''iO,K;i) 



j=0 



< ln + 



j=0 



< In + Y^NnJ ulj \\+M4 deulj I 

< / 
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where Kr, is a constant which only depends on C2, Km and K]\f. If we pick H larger than we 
get 

Finally, we estimate r/^ using point 3. of Lemma lOl and point i. of Lemma [6^ 
iK.+ill < ir4/nC22"r„,^(2")|li + M„(c22«r,,^(2«)|liy 

< ^ + 



< 



16 2» 2*16 2* 

Wn 



2i+i- 

We used here that K > K4C2 and K > KmCI- 

End of the first part. We apply the operation "solving the uq equation" until the size of 
ll^nill becomes smaller than '^"+^ ; more precisely, when ~ is smaller than "''^^^ . At this time 
in we stop and we pass to the second part of the stage. 

7.4.2 Second part: the ni equation 

This part consists in solving the equation for ui only once. More precisely, we search for a 
coordinates change of the form h\^{6,z) = {6,Ui^{9)z) such that 

^;||f^{.?:rw + A} = A^. (68) 

Notice that if we put e"""^^^ = Ui^{6), equation (|68p reduces to solving 

vn9)-v^{9 + a)=logXi-log{p'l;i"{e) + X). (69) 

However, due to technical reasons we will just deal with the following truncate equation 

VtiO) - + a)= log - { log {pli^ie) + X)}\^„. (70) 

The complex number is the unique one for which a continuous solution to this equation exists 
(see Section 15. 3p . More precisely, it corresponds to the unique complex number for which the 
mean of the r.h.s. vanishes, that is, 



log Xn 



We have 



j^^ { log (p^;;"(0) + x)}\^^de = log A + J^^ log + i)de. (7i) 
j^\og{Pi^ + i)de\<c,\\pr"\\ 



for a constant C3. To estimate the solution of (ITOp we need an estimate on the size of the r.h.s. 
For this, we rewrite this expression as 

log A*. - log (p^f (^) + A) = log + l)d0 - log (pj- + 1 
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and we obtain 

||iogA^-iog {p'l:t-{e) + x)\\<c4Pi:n\ 

for a constant C4. We have 

ll^^ni < C2C^2^r^{2^)\\plt\\ (72) 
\\devn < C2C^rr^{2^)\\pli"\\ (73) 

since we choose the zero mean solution. Letting = exp(u^^) and conjugating by /i^ = 
(6',u"j(0)z) we get the new functions 

The last one is defined on a disc D(0, r"^^''^(0)) for the z variable, with 

KM 



= (76) 



The sequence {i?"'*} verifies 



for every in the strip Bgn . From now on, the domain for the z variable is the disc -D(0, K 



•n+l,0\ 



The untouched rest in equation (I70p gives raise to a new function pi defined by 
Estimates at the second part of the stage. By (172]) and (f6T]) we have 

- 111 < C3C2C42"T,(2")K4/n. (79) 

Therefore, there exists a constant -R'e such that 

\\ul,\\ < e^«2"r.(2");„ (80) 
||«,,)-i < e^«2"r.(2");„_ (8^) 

Using the equality dgu^ = dgv^Ui and ([73]) we obtain 

WdeultW < e^«2"r„(2")/„^g4np^(2n);^ (82) 
< 3i^64"r,(2")Z„ (83) 

provided that H is larger than Kq. We remark that the estimate (|79|) says that the topological 
degree of ui is zero. We may now proceed to conclude the definition of the sequence {i?"'*} and 
to prove the lower bound 3/8 for it. 
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The sequence Assume that R'^'°° is aheady defined. We put 

nn,oo 

nn+1,0 _ ^ foA\ 

^ g«2"r^{2")/„ • 

With this definition, inequahty (j77p fohows from (j76p and (|8ip . The sequence i?"'* also satisfies 
the monotonicity properties announced in (j36p . Finahy, to check that i?"'" (and therefore each 
ii""'*) is bounded from below by |, we first note that 

n 
j=0 



Each of the differences above can be estimated by 



2C722^r,,/3(2^>,- 



g«2JTa(2J)«^ 

< K2^r,(2^)/j + 2C22Jr„,;5(2^>,. 
Thus, from point 2. of Lemma 16. II and point 2. of Lemma 16.21 

(iJ^'O - < ^ H2"r«(2")/„ + 2 ^ C22"r„,;3(2")«;„ < - 

n n 



Relation (j75[) gives 

o2 n+1.0//i \ r,2 n.in/n n fn\ \ 

Therefore, according to the definition below of the sequence {M„} we get the desired estimate 
for d^pt, namely, 

= Mn+i. (85) 

The sequence {M„}. We first let Mn* = M, and assuming that Mn is defined, we let 

M„+i = M„e3«2"r.(2'^)^., (86) 

According to this definition, one easily checks the existence of an upper bound Km for the 
sequence M„ ■ 

In order to compute the oscillation dgdzP^^^'^ we first note that 

dedzPt {0,z) = ^ dzPt' {e,u^t{e)z) 



dedzPr"{0,nUe)z) 



22 



which imphes 

osc(a,9,pr^'°) < 4||5e<J||<J||KJ-i2i^M + A^n,J|<tllllK,)-i 

+2\\{ul,)-^\\ulMM\\deult\\. (87) 

At this point we can conclude the definition of the sequence {Nn^i} and show the existence of 
an upper bound Kj\f for it. 

The sequence {Nn^i}. Assuming that Nn^^o is aheady defined, we put 

A^„+i,o = e2^«2"r'^(2")/„^^^^ ^ 144KMi^64"r,(2")/„. (88) 
This definition together with inequahty (I87|) and other estimates in this Section give us 

osc(ae9,/jr^^'°) < Nn+1,0. (89) 
By replacing the value of Nn^oo (see ((671) ) in the definition of A^^^+i^o we obtain 

iVn+1,0 < e2^«2"r.(2")«„^^^^ ^ 8KMC24"r„,^(2")u;„ + 144i^Mi^64"r«(2")Z, 

< e2-^«2"r.(2")«„^^^^ ^ i^74"r„(2")/„ 

< e2^« ^2"r„(2")«„ j^2iV + 5] 4"r„(2")/„) . 

These last inequalities hold for a well defined constant Kj. We see that (|9Up and the properties 
of the sequence {In} imply an upper bound Kj\f for the sub-sequence Nnfl- Finally, due to the 
following monotonicity relations, 

2N = iV„.,o < Nn*,l < ■ ■ ■ < A^n*,oo < A^n*+1,0 < • • • 

this implies the existence of an upper bound for the whole sequence {Nn^i} ■ 

We conclude this second part of the stage by writing out the functions whose expressions 
involve some truncate functions and for which we need to allow a loss in the width of the strip 
in order to get good estimates. First of all, relation (j78p implies 

dep'li''' = A* (e^-(^"^'^+^)}l^)a,({log(,-^ + A)}| J 
Since the truncation operator commutes with the derivation with respect to 9, we have 

Finally, the expression for pQ is 

= ^^ip;r^) = ^lyjT^ [<^^ + E ■ (91) 
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7.4.3 Third part: loss on the width of the strip 

In order to estimate the size of the truncate functions appearing in PU]) and ([HT]) . we introduce 
a loss in the width of the strip of size max(d^, d^). The new domain for the 9 variable is (5""''^. 
We can estimate Po^^'" by 



II n+l,0|| ^ o / , \ - ^ Wne 

IK* II < 

< Wn+l 

provided that H is larger than 2Ci. Note that here we used the estimate ||('Uit)~^|| < 2. This 
estimate is not sharp but it is enough for our purposes. To estimate the expression in (f90l) we 
first notice that, since ||/Oxt"ll is uniformly small, there exists a constant C5 such that 



C5||5ePi'M|e- 



-27r2"d 



n j1 



< 



< In+i. (92) 

This last inequality holds thanks to the definition of the sequence {ln\n>o (see Lemma [6. ip and 
the fact that H is greater than 2C5 . In this way, we get the family of inequalities (I38p . . . (I4ip for 
the index n + 1. This finishes the third part. 

7.5 Estimates with respect to the parameter 

Before dealing with the fourth and last part of the stage, we treat again the first three parts in 
order to get some estimates on the parameter t. These estimates will be essential in the fourth 
part, where we will reduce the parameter space in an important way. However, in this Section 
we will also reduce the parameter space in order to get some useful Cauchy estimates. At the 
beginning of the stage n we have the following bounds: 

osc(ata,p;^'°) < r„,o (93) 

atp^fd^-Ao < Sn. (94) 

Ti 

The sequence of real numbers {Tn,i\ n>n* ,&i\j{oo} verifies 

Tn+ifl < i^Tn2'^r,(2") (95) 
for a constant Kt- The sequence of real numbers {sn}n>n* verifies 



< < Sn*+2 < ■ ■ < -j-. (96) 



1000 " " 100 
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Estimates at the first part of the stage. We will show inductively on i > that 

osc{dtd,p''/) < Tn,i. (97) 

Suppose that this estimate holds for < j < i + 1. We start by estimating dtu^'l- For this we 
use a Cauchy estimate, that is, we lose a small quantity in the radius of the disc D{tn,Pn) (the 
parameter space) in order to get an estimate for dtu^'l- More precisely, we construct a sequence 
of radius {p ~ } jeNuoo defined by 



Pn- = Pn~-^n 

I z — 1 

Pn- = Pn- 3ln- 



2\i 

3. 



This says that at each step we lose a very small quantity of radius of size ^n(^|j (we use this 
loss in the estimate (|98|) that follows). We note that at the limit, the total loss has size 3ln- 
From (I62p the Cauchy estimate gives 

iia«siii<g(j)'. m 

We can now estimate the oscillation of dtdzP by 

dr\ n.i-\-l o n.i / /\ n.i , \ , o2 n.i / /\ n.i , \o n,i 

td.pt = dtd.pt {9, Uq'j + z) + d,pt {0, n^'j + z)dtUf^'^t 

tOzPt {0,UQ^)-d^Pt {e,UQ^)dtUQ 

OSc{dtd,pf+') < Tn,i + 2Mn\\dtU^i\\. 

This enables us to define the sequence {Tn^i}i>o by letting T„*^o = 2T and defining recursively 



N2" V4 
We define T„^oo by 

T„,^ = T„„o + ^^. (100) 
Notice that Tnfl is an upper bound for the Tn^i, and the following inequality holds 

osc{dtdzP^''~^^) <Tn,i+i. 

Estimates at the second part of the stage. First, we estimate the size of dfU^ = uldtVn- 
By taking derivatives with respect to t in both sides of ([70]) we obtain 



d,vm-dtv^{e + a) = -{^2^\l . (101) 



Pl,t +^ 

Recall that the notation | f " means that we only keep the Fourier series terms with orders between 
1 and 2". The best estimate that we can get for the right hand side function (in his original 
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untruncated form) is of the order of a constant. More precisely, by introducing a loss in the 
radius of the disc D{tn,p„-) of size Ifi , W6 get a Caucliy estimate 



tPit 



< l2Ki. 

Pl,t + ^ 

The new parameter space disc is D{tn,Pn-)j where 

Pn~ — — P''^ 

Then we have 

WdtultW < C224K42"r„(2"). 
We may also estimate the oscillation of dtdzPt^^'^ by 



(102) 



(103) 



dtdzpl 



+ 



ui,{e + a) 
+dlpr%e,uim^)^dtui,{e)} 



o.c(9i9,pr^'°) < 4||9i<,J|||<,J||(<,)-if M„ + ||<,J||(<,,)-i||r„,,„ 
+2\\ulMWlt)-'\\Mn\\dtul^\\. 



(104) 



We introduce now the definition of the sequence {Tnfl}n>n* and show the claimed upper bound 
(ESI). 



The sequence {T„^j}. Suppose Tn^oo is defined. We put 

Tn+1,0 = i^T2"rQ,(2") + 2Tn 



(105) 



where the constant Kt is given by Kt = 32KmC2'^4:K4^. By looking at the inequality ()104p we 
obtain 

If we replace the value of Tn^oo (see (jlOOp ) in the definition of Tn+ifl we get 

Kt 

Tn+1,0 < i^T2"r,(2") + 2Tn,0 + 
This gives the desired estimate (in fact, we get a sharper estimate) ■ 

Estimates at the third part of the stage. We can estimate the distance between Jji dtpi^t 
and Aq and show that this distance is small. In Section 17.61 this will allow us to show that 
jrpdtpi^t{Q)d6 grows almost like a linear map. In order to simplify the notation, we introduce 
the constants i^g. Kg in such a way that 



dtij^dzp'l'^e^uim^ < T^,oo2C22-T^^p{2-)wn + ^ 



< 



(106) 
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We will also need some estimates on the size of (9t(A^). For this, it will be very useful to have a 
more developped expression for it. From (|7ip we have 



xidtiiogxi) 



A 

X* 

X 



By differentiating (jTHjl we get 



■de 



+ A J 



The difference between the integral of dt(\[^'^ and Aq may be written as the sum Ii + I2 + /s 
of three terms corresponding to the expressions in the three lines below: 



/ 



yt g{log{p^''"+A)}|2^ 

A 



dQ 



QtP\,t 

n,in 



d9] - An 



+ 



A 



de 



dtiVo-'dzP'^'iO,^:})) 



de 



1 + 



Pl'.t 



+ / A*e^'°s(prr+A)}l2^/^!f!l:^\ 



dn,t, 
tPu 

Notice that the term Ii itself equals the sum of the following three terms 



{log(p^-'"+A)}|2^ 

A 



de 



tPi 



n,0 



1 



+ 



A 



J log{^+l)+{log{^^+l)}|2 



^0- 



Hence, there exists a constant Kg such that 

llllll < e^^^"Sn + Ksln- 



(107) 



The terms A^,p"'^*" are bounded from above by a universal constant. This together with ()106p 
implies that there exists a constant Cq such that 



II/2II < 



H2" 



(108) 



The term I3 equals the integral of a universally bounded term times the truncate of a function 
whose size is also bounded by a type K constant. Notice that these bounds exist even before 
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the introduction of any loss on the width of the strip, see ()102p ). When introducing the loss on 
the width, we obtain a constant Kiq for which the following bound holds: 

By considering these three bounds for the integrals Ii,/2 and Is, we get the following estimate 
at the third part of the stage: 

X. ftft-f '-^" - ^»|| + + ^ + iJSr- ("0) 

In the next paragraph we will give an explicit definition of the sequence {sn} and also prove its 
claimed properties. 



The sequence We define this sequence by recurrence by letting 



1000 



(111) 



= ^ne-^'"+i^A + ^ + ^. (112) 
Relation (jllOp implies (tautologically) that 

/ 5tp^f'°de-Ao <Sn+l. (113) 
By iterating the definition we get the estimate 

OO OO OO 1^ 





OO 



2C6K9 , 



10 



u 

If is large enough then this gives (j96p as desired ■ 

7.6 Fourth part, reduction on the parameter space 

Now we arrive to a different stage of our process. Indeed, the estimates here do not arise 
neither as a consequence of the remainders of a linearized equation (quadratic residue) nor as 
the result of some Cauchy estimates associated to some loss in the strip. In this Section, we get 
our estimates as a consequence of the transversally properties of the family {Ft}. Let's see the 
necessity of such a hypothesis. From relation ()78p . at this moment we have a constant C7 such 
that 

\\p'l^'''\\<Cr\\pli"\\<CrK,ln. (114) 

In other words, the size of p"^^'*^ is still of order /„. This fact is not a surprise, since we have 
made operations which deal only with the size of the derivative dgpi. We will use the Affirmation 
IS.ll to control pi. Thus, we need to control the size of the complex number /^i pi^t{G)d6. To do 
this we will strongly use the transversality property for the family {Ft}. We will find a parameter 
tn+i such that fji Pi~l^^f^ = 0. Then, we will exploit the fact that the size of f^i Pi~l^^'^ is very 
small if we are placed very near this special parameter, due to the continuity of the functions. 
More precisely, as in ()37p we will reduce the parameter space D{tn,Pn-) by finding a new disc 
D{tn+i,Pn+i) C D{tn,Pn~) with i„+i being a simple zero of f^.^ Pi~l^^'^{6)d6 and Pn+i small 
enough. 
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7.6.1 An application of Rouche Lemma 

We will use the classical Rouche Lemma to find tn+i- To do this, we need to show that the 
hypothesis of the Lemma are satisfied on the boundary of some disc around We develop 

= j (e^Ti log de+{iog(pJ;^ {e)+mw - 1^ dO. 

We then define g+{t) as being the value of this last expression. We compare this function with 
the function g(t) = Jji — dO on the boundary of some disc D(tn,r): 

\g+{t)-g{t)\ < ^^||^^log(^l^ + l)d^-5(i) + {logK,f(^) + A)}|2^ 

where Z = \ j^^ log (^^^^ + l)de + { log{p1'l"{e) + A)}|27r|. Then there exists a constant 
Kii such that Z < Kuln < 1, provided that is large enough. Thus, using the Taylor series 



expansion of log I h 1 1 , we have 



\g+it)-git)\ < 



dO 



+XJ{'"«5(''m"w+^)} 

+ Xj{log(pr,fW + A)} 



2" 



de 



2" 



2 V A 

. + Z^e^ de 



de 



where the absolute value of oo{e,t) lies between and 1. One can easily check the following 
estimates: 



/ \ Y.d,pr{e,ull{e)) 

2 \ \ ) 

J{iog(pi:rw+A)} 



de < 



2C2KMln 



de < ^ 

2 

at/ < < 



di 



Therefore, if M: is large enough, for every r < p„- and every t G D{tn, r) we have 

\g+{t)-g{t)\ <3lnL-\ 



(115) 

(116) 
(117) 

(118) 



The transvcrsality condition on the derivative dt /yPi,t(^)(i^ implies that g{t) grows almost as 
a linear application, namely, 



I ^ I — 



Ao(t - 
A 



g{t)+gitn) 
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This implies that the inequahty 



99 99 

\g(t)\ > (L^^ - Sn) \t - tJ > \t - tJ = r 

holds on the boundary of a disc D{tn,r). Therefore, if the radius r is at least ^3/^, then the 
hypothesis in the Rouche Lemma is satisfied. So, we get a simple zero, tn+i, for inside 
D{tn,r). Now we need to show that if we choose r as being this minimal value, then the new 
zero tn+i lies effectively inside the parameter space D{tn,Pn-)- To do this, it is enough to see 
that 

Pn- - > 92/„ > 0. (119) 

Moreover, we need to show that we can start with the stage n + 1 having a parameter space disc 
of radius Pn+i- For this we notice that 

Pn+l = 100/rt+i 

ioor,,^(2-)r„(2") 



^ra,/3(2"+i)r„(2"+i)2e2'^+i . 
< 25Z„. (120) 

From pi9|) and p20|) we conclude that the disc I?(t„+i,p„+i) is contained in Z)(t„,p„-). Thus, 
the new parameter space is the disc D{tn+i,Pn+i)- 

In this way we get the whole necessary estimates in order to relaunch our process at the stage 
n + 1, and this finishes the stage n. 

7.7 Convergence of the method 

As we have seen in the previous Sections, we can iterate the algorithm infinitely times for the 
fhd corresponding to the parameter t = r\n>n*D{tn,Pn)- (Of course, t is the parameter claimed 
in the Theorem 14.11 ) In what follows, we will omit this parameter in all of the notation, and 
we will deal with F = Fj. We need to show that the successive compositions of the conjugacies 
constructed in the course of our process gives raise to an analytical change of coordinates. The 
decreasing rates for the functions and all of the estimates arising in our process are sufficient for 
this. We let 

OO J- 

<5- = 5*- j;max(dO,4)>-, (121) 

n* 

and we denote 

Hni9, z) = hnfi o hn,l o • • • o /l„,i„_i o /l„(6', z) = (^6, Uq{6) + Ui{0)z^ , 

where we write Uq{6) = Yli^^ ^o'*(^)- The successive compositions of the conjugacies up to the 
stage n is 

Hr,* O Hn'+i 0...0 Hnie, z) = {6, uf + uf (u^'+l + uf+\. . . {u^ + U^^z) ...))) 

/ n-1 \ / n 
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We need to show that the constant term and the coefficient of z above, converge uniformly 
on Q when n goes to infinity. For this, it is enough to see that the product u\{6) converges 
uniformly and the series ^o(^) converges absolutely and uniformly. For the product we have 

\j=n* I j=n* 

The above series converges absolutely and uniformly by the estimate (I72p and point 2. of the 
Lemma |6. 11 The series of uq is 

n n 
j=n* 3=1^* i=0 

which allows to obtain the estimate 

j=n* j=n* i=0 j=n* 

Thus, the following limit exists 

H{9, z) = lim Hn* o Hn*+i o ■■■ o Hn{0, z). 

n—>oo 

We write this limit in the form 

H{9,z) = {e,uo{9) + ui{9)z), 

where the functions uq,ui are analytical on the strip Bs°°. The function ui has zero topological 
degree since the topological degree is a homomorphism and each function u\ has zero degree. 
The conjugacy F hy H takes the form 

H-^oFo H{9, z) = {9 + a, Xz + p°°{9, z)). 

Indeed, the sequences ln,Wn which control the size of Pq'^,Pi'^ go to when n goes to infinity. 
The function vanishes up to the order 2 at z = and it is defined on the disc D{0,R°°), 
with I < R°°. The invariant curve is uo{9), and due to ()122p it verifies 

1 °° 1 

WuoWbs^ < - ^ — ^. (123) 

j=n* 

7.8 A modification at the first stage and tfie size of e 

Let n > n* be a natural number and e be a real positive number such that 

— ln+1 <£< —h- (124) 

Let {Ft\t&T, be an analytical family verifying the hypothesis of our Main Theorem and 

„ „ „ „ 24e , , 

We start the algorithm at the stage n. We notice that in order to pass to the stage n + 1 we 
need to show the following facts at the end of the stage n: 
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In the first part. The inequalities (j59l) . (l60|) must be 

\\depf;\\ < 2(^). (127) 
In the second part. Inequality (|79|) must be 

24e 

Now the first two inequalities follow by repeating adequately the computations of the first part 
while the third one follows immediately. 



ult - 1|| < C3C2C42*'r„(2")i^4 — • (128) 



In this way, we reach the third part and we obtain successfully the corresponding bounds for 
/fi+i and Wn+i when introducing the loss in the strip. 

In the estimates with respect to the parameter. In order to get ([981) it is enough at 

each step to alow a loss of size ^ for the radius of the parameter space disc. Moreover, 

to obtain p(J2|) it is enough to alow a loss of size ^ . At the fourth part we can see that (jllSp , 
()116p . (IllTj) . are bounded from above by 

^6 ' 2 ' 4 « ^^^^^ 

respectively. If H is large enough then (|118p is bounded from above by 3L~^(^^^^. This allows 
us to finish the stage with an appropriate loss of width for the strip. 

The above considerations show that in order to start with the stage n we just need a param- 
eter space disc of radius 

24e , , 

Pn = 100—. (130) 

For such a disc, we also verify that the parameter t proving the existence of the invariant curve 
is located at a distance e times a constant from the origin in the complex plane, as desired. If e 
goes to 0, then the level n of the first stage grows and (I123p allows to conclude that the size of 
the invariant curve tends to 0. 

Actually, we have shown that Theorem 14.11 holds under the stronger hypothesis \\po\\ < Wn* 
(or more precisely under hypothesis (I125p ). In the next Section we will see that we only need 
to impose the hypothesis ||po|| ^ £• We will see that even under this weaker hypothesis we can 
reduce the problem to the stronger hypothesis used up so far. Indeed, we will describe a previous 
preparative process showing this reduction. 

8 A previous preparative process and the end of the proof 
8.1 A initial dicrease in the size of po 

We know that solving the uq equation gives rise to a diminution on the size of po. Thus, we will 
solve many times this equation to get the adequate size for po in order to start the algorithm. 
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Although each time we solve the equation we slightly increase the size of pi, the total amount of 
the increase is well behaved (and may be effectively controlled). We will omit many details when 
manipulating these estimates, the corresponding justifications being contained in the previous 
Sections. 

We fix a constant by using the data (^L,M,2T + {2000L)~^,d, as described at the 

begining of the Section EH This choice gives us our sequences {In} , {wn} , {d^} , {dn} and a 
natural number n* representing the lowest starting stage for the algorithm. Suppose that there 
exists n > n* and e > such that 

-^<eM <^. (131) 

Assume also that the hypothesis of Theorem 14.11 hold and that for every parameter t in the disc 
D{0,Kji£) we have 

\\po,t\\ < £ , \\pi,t\\ < e, (132) 

where Kji is a constant which will be defined in a sequel. The value of the radius for the 
parameter space disc corresponds exactly to the sum of what we need to start this previous 
preparative process and the necessary radius to start with the algorithm at the stage n. This 
previous preparative process also has stages (indexed by the natural numbers between n* and 
n — 1). Each one of these stages contains two parts. The first part of the stage n is at some 
moment divided into steps. Each step consists in solving the equation of uq once, with truncation 
up to the order 2". The second part consists in introducing a loss for the width of the strip of 
size d^. For notation reasons we will use the sequence of real numbers 

Wn = C2M2"r„,^(2"). (133) 

We present the first stage n* separately. The reason for doing this is twofold. On the one hand, 
at this stage there the estimates are slightly different. On the other hand, this stage allows us 
to understand the manipulations and estimates that arise in the remaining stages. As usual (in 
this work), we denote the starting functions as p^\, pi\, pf* . We introduce an index to indicate 
the steps. We define the sequence of functions r/^, ^ by 

vi*,o = Po/ (134) 
Vn*,i+1 = Pl,t ^o,t +Pt (^^o.t )• (135) 
The hypothesis of Theorem 14.11 sav that ||r/^* q|| < e. We will show inductively that 

II t II ^ K/n*,0 , £ rioc\ 
\\Vn*,i\\ ^ -^^<Ji (136) 

^ \\Pi7\\ + 2|l^n*,ol|W^n* < e(l + 2Wn') < e{wn*)-\ (137) 

To do this, assume that these inequalities hold for < j < i + Then we can estimate the new 
V by 

<e(l + 2Ty„,)hi.,JI^+M(||r?*., " 



M 



Since e{l + 2Wn*) < £{wn*) ^, ||??^. ill ^ ^ n> n* we have 



Wr,* dW„* 1 , W'^, W'^, 1 
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This implies that ||r/^* < ||r/„*^j||2 ^. The size of H/o" ^'^^"'^11 is bounded from above by 

i 

\\pf/^'\\ < WptAl + Wn'Y.W'In^J < \\Pl7\\+'^Wn4vi*,0\\<e{l + 2Wn*). 



We continue the steps until to the moment when the size of rj is smaller than or more 

precisely when ~ < ^^^'^^ (this subtle difference will be very important in Section [9]). Let in* 
be such a moment. In the second part of the stage we introduce a loss in the width of the strip 
of size . The final size of po is thus bounded by 



^0,t^'^\\ ^ ll^n*,j„. II + 2||??n*,oll 



, e 



-27r2" rfO. 



d2,* 



< 



ewn*+i 



8.1.1 An arbitrary stage n* < n < n — 1 

At the begin of the stage we have 

KfW S ^ (138) 

n,0|i / , oTTA , '2Wn*+lWn*+l , , 2Wn-lWn-l 



IK?II < .[l^2W„.^ + ,139) 

< e{wn*y^. (140) 
For the usual sequence {r/^ j} we will show inductively that 

IKJI < ^ (141) 
< e(l + 2Wn* + --- + ^^^)<e{w„,r\ (142) 
So, suppose these inequalities hold for < j < i + 1. The size of r/ is bounded from above by 

^(-n.)-^^ll<.ll+M(||<,||^)'. 

Moreover, we have 

eiWr,*) < tn < — 

\ n ) ^ n24M 4 

That implies ||f?^ j_|_il| < ||?7^ jp^-*^. The size of Pi'l~^^ is bounded by 

WpT'w < ii<ii + ^"Eii<J 



^ II n,0|| , 2WneWn 
< "K.ll + ^^ 
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Wc continue the steps up to the first moment where the size of r] is smaller than 2Hw^I ' °^ more 
preeisely when ^^"^'"^^ < ^jj^- Let i„ be sueh moment. At the second part of the stage we 
introduce a loss in the width of the strip of size dj^. The final size of po is bounded by 

-27r2"(i0 

IKfll < ll<^JI+2||<ol|— ^ 

^ £Wn+l _^ ^_£Wn £Wn+l 



< 



£Wn+l 



Hence, at the end of the (n — l)-th stage we obtain 

II n II < ^ 6e{Wn')-' , . 

"1^°'*"" - < H4-r„,^(2-) ^^^^^ 

WpltW < ^KO"'- (144) 

These functions are defined on the strip S — Yln''^ "^^^ function Pt{9,-) is defined, for the z 
variable, on a disc of radius 

n-l in 
n=n* i=0 

If we start the algorithm at the stage n with this parameter space disc, we will reach the end 
of the process with a radius even larger than 3/8. The size of does not change with the 
operation of resolution of the equation of uq. 

8.1.2 Estimates with respect to the parameter 

We will need to estimate the values of the following two sums: 

n-l in n-1 in 

EE Kill ' EEii«:^ii- 

n=n* 1=0 n=n* i=0 

The first one is bounded by 

n=n* 1=0 n>n* 



2Wn*e ( 1 \ ^ 4W„* 

n>n* " 



By introducing a loss of size 96000LW„*£ on the radius of the parameter space disc, we obtain 
a Cauchy estimate 



n-l in 



n=n* i=0 

Now, we can estimate the final oscillation of dtdzPt by 

n-l in 1 

osc{dtd,pf) <2T + 2MY,Y. IK:* II ^ 2r + 2000^. (147) 

n=n* i=0 
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We will need the following estimates 



n—l i 



n=n* i=0 



< 



< 



1 N 4Wn* 1 

2T -I I — e -I 

2L) M 4000L 

1 



2000L 



where the last inequality follows since the size of ^^ is very large. We define the complex number 
Ao by 

(148) 



t=to 



8.2 A simple zero for J^^ pi t 

There is only one remaining ingredient we miss in order to apply the algorithm to our family, 
namely, the existence of a simple zero for jr^i p^^ at the center of the parameter space disc. We 
will find this zero by using Rouche Lemma. Thus, we need to compare J-jpi p^^ with the linear 
part of /rj,! p"^ J at to on the boundary of a disc of radius R and centered at Iq : 



Plt-dti / Pit ^ ^ it -to) 

Jl \ Jjl / t=to 



< 



+ 



Plt-dti / Pi,t] ^ ^ it -to) 

Jl \ Jjl / t=to 



Pit - / Pit 



T 



< e + -R' + 4Wn*e. 



We have the linear estimate 



dt[ I Pi,t 



over this boundary. We choose the radius 



t=to 



it - to) 



Rr 



L-1 - v/^-2_2re(l + 4Ty„0 



(149) 



which corresponds to the smallest root of the equation between the two precedent bounds. In 
this way we assure the existence of a simple zero t* of Jji p^^ inside Dito,Re)- We easily verify 
that i?£ is a positive real number. Furthermore, we also verify that 



(150) 



because a; > 1 — \/l — x for every small x. 



In order to start the iterative part of our process (the algorithm described in Sections 
17. 4[ . . 17. 6p we need a parameter space disc of radius 100 ^^^^"'^*^ — and centered at t*. All what 
we have seen allows us to define the constant Kr as being precisely what is needed. That is, to 
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find t*, and to lose a little bit of radius to obtain the estimates with respect to the parameter 
in the previous preparative process. Thus we have to choose 



Kr = 2L(1 + iWn') + 96000LW„* + 



2400 (w„ 



(151) 



To finish this previous preparative process and to start with the algorithm, it is enough to notice 
that the following estimate 



^o-dt[ I Pit 



< 



^o-dt[ I Pit 

n—l i 



+ 



n=n* i=0 



< Te(2L(l + 4iy„.)+ ^''°»'r''' U^'' 



2000 



< 



L- 



1000 



holds inside the parameter space disc centered at to with radious e( 2L{l + 4Wn*) + ^^'^Q^^"* ) — 



Notice that this disc contains the disc centered at t* with radius e 
The invariant curve is 

n—l in h n 



2400(t»„,)-i 



n=7i* j=0 



(152) 



]=n 



It has zero degree because all the coordinate changes performed in this previous preparative 
process have zero degree ■ 



9 A parametrized version of the persistence of the invariant 
curve Theorem for the fibered holomorphic dynamics 

Let A G C be an open set. We consider an application from A to the set of analytical 1 complex 
parameter families of fibered holomorphic dynamics, 

sEA^{F/}ies. 
We say that this application is analytic if 

{s,t,e,z)^Ft{e,z) 

is an analytical function. Let s be in A. Suppose that the family {-F/}tes verify the hypothesis 
of Theorem 14. 1[ Hence, for every t in D{tQ, Kjie) C S one has 



L > 



9t[ /ti pU 



t \t=to 



< £ 



< e 



\dip!\\<M , \\dtd,pt\\ + \\diplt\\<T 
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for some e € (0,e*], where e* is given by Theorem 14.11 Then, there exists a neighborhood V{s) 
of s in A such that the following holds: every family {Ff}t^Y,, with s in V{s) verifies the above 
hypothesis. Theorem 14.11 gives us an application from V{s) to D{to, Kre) 



and an application from V{s) to the set of analytical curves from the circle to C 

s ! — >u'' :Bs^C. 

The curve is invariant by the fibered holomorphic dynamics , has zero degree and its 
fibered rotation number equals /3. The goal of this Section is to show the 

Theorem 9.1 The applications s ^ tg, s ^ defined above are holomorphic functions. 

The proof is obtained by making a more carefully study of the method (previous process and 
iterative algorithm) , and more precisely by the ability of performing this method in an uniform 
way for the whole set of families {F/j^gs, s £ V{s)- Let n be the unique natural number 
verifying 

We can apply the previous process to every family {Ff}t£-s, s in V{s). This gives us the 
functions {(mq ^ )^}(n*<n<n;0<j<i„) defined for every t belonging to D{to, Kre) (here the index 
are uniform for every parameter s in V{s)); this also gives us the functions {pQtY, (PitY, {PtY • 
These functions have an appropriate size to start with the iterative algorithm at the stage n. 
Furthermore, these functions are defined for every t in D{tfi,pn). The parameter is the only 
simple zero for Jji{p^ ^.ydO given by the last operation in the previous procedure. Each stage 

n of the algorithm (n > n) produce at some moment the functions {{uQ'ly}o<i<i„, (uit)'^ (with 
uniform index). These are defined for every t in D{tf^,pn), where t^ is the only simple zero for 
IjiiPitydO. The algorithm also produces the functions {Po'l^Y, (Pit^Y, iPt~^^Yj which have an 
appropriate size in order to start the stage n + 1. These functions are defined for every t in 
-D(i^+i,Pn+i). 

Note that from the stage n and so on, the parameter space disc is the only thing that is not 
uniform (the center is not uniform). We will see that all the above functions are holomorphic 
functions with respect to the variables {9, z, s,t). 

9.1 The apphcation tg 

Note first that if one consider only one family {Ft}, the parameter t for which one finds the 
invariant curve equals the only element in the intersection Un>nD{tn,Pn)- As p„ ^ when 
n — > oo, one deduces that t equals the limit of the sequence tn of the centers of the discs. This 
suggests a strategy to show that tg is holomorphic. We need to show that each application 
s I— > is holomorphic. In this case, the (uniform) limit tg is also holomorphic, since the radius 
Pn are uniform. 

To show that s i— > is holomorphic we use a following well known fact: Let {gg : B — > C}sgA 
be a family of holomorphic functions depending analytically on a parameter s and verifying that 
each function gg has a unique zero in D which is simple. Then this zero is a holomorphic function 
on the parameter s (this follows directly from the classical Implicit Function Theorem). In the 
next Sections we show that the function Jji(p^ ^Yd9 is holomorphic on s and t. 
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9.2 The functions (u^'lY, for n* < n < n 

For every s in V{s) the notation {uq'IY has an uniform sense with respect to the index, provided 
that n* < n < n. This denotes the solution of the equation of uq for the 2"-order truncate of the 
function (r/^ J*. Then {uq'IY is a trigonometric polynomial of degree 2". Its coefficients depend 
linearly on the Fourier series coefficients of (ry^ J**. Since (ry^* qY is holomorphic with respect to 

s and t, the function (uq^''^)* is also holomorphic with respect to s and t. 

We denote by < the (direct) lexicographical order on N x N. Notice that an index {m,j) 
appears in the algorithm before another index (n,z) if and only if (m, j)<{n,i). Let {n,i) be 
an index. If for every index (m, j)<{n,i) we have that {i]m jY holomorphic with respect to s 
and t, then {uq'^Y is holomorphic with respect to s and t. By an inductive reasoning we obtain 

that all of the functions {(^^o j )'*}{n*<n<n;0<i<j„); ^ '^^11 functions (Pq {PitY ■• {PtY ^-^^^ 

holomorphic with respect to s and t. In particular, the function J,jj,i (p" is, and consequently 
t* = i^, is holomorphic with respect to s. 

9.3 An arbitrary stage n > n 

Suppose that for every index (m,j)<{n, 0) all of the functions (PotY^ (PitY^ {PtY ■> {i'^^tY} ^'^^ 
holomorphic on s and t . Suppose also that for every m < n the centers of the parameter space 
discs are holomorphic functions on s. For every s in V{s), when solving the first equation of uq 
in the stage (the step i = 0), the trigonometrical polynomial solution {u^fY S'lso depends in an 

holomorphic way on s and t. As a consequence we have that the functions {p^'^Y^ (Pi't)^^ (p"'^)^ 
are holomorphic. By repeating the argument we can get the same conclusion for every function 
W)o<i<i„' (^i,*)''- This implies that the same holds for the functions (/>o,^^)^ {pi^^Y, (.p7^^Y- 
In particular, the functions 




(153) 



which are holomorphic on t £ D{t'^,pn), are also holomorphic on s. Then there exists a holo- 
morphic application s ^ t^^^ such that 

l^yi,t^^J^de = 0. (154) 

In this way, we show inductively on n that the center applications s tf^ are holomorphic. Thus, 
the limit application s i— > t ^ also corresponds to a holomorphic function on s. The invariant 
curve Us is written as (see I152|) 

n— 1 i„ h n 

n=n* i=0 n=n j=n 

The first term corresponds to the previous preparative process. This is a finite sum of holomor- 
phic functions on s. The functions arising from the iterative algorithm are also holomorphic on 
s. We may then conclude that the application s i-^ is holomorphic. This completes the proof 
of the Theorem 19.11 ■ 
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